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1. Introduction 


Future gravitational-wave astronomy requires gravitational-wave detec¬ 
tors of fantastic high sensitivity. Two designs of detectors are on the beam of 
realisation: bar detectors and laser-beam detectors. Both are promising to be¬ 
come enough sensitive for the detection and even for the detailed measurement 
of cosmic gravitational waves. Whereas the beam detectors will reach rather 
high sensitivity the bar detectors will be quite cheap. The sensitivity of the lat¬ 
ter should be sufficient to detect gravitational waves from coalescing binaries in 
the VIRGO cluster and from type-II supernovae in our galaxy. Unfortunately, 
both of them happen not very often (about one per 30 years). However, there 
might exist other sources in our galaxy, like slightly deformed rotating neutron 
stars, in particular young asymmetric pulsars, or even completely unexpected 
sources which make gravitational-wave surveys in our galaxy, or other galaxies 
of the Local Group, very desirable. The centers of galaxies are very difficult 
to investigate with present-day telescopes so that our knowledge of the centers 
is very crude. Gravitational-wave detection might reveal interesting objects 
there. 

Purpose of this paper is to work out in details a general formalism for 
treating the interaction of a gravitational wave with an electromagnetic (ohmic 
and superconducting) device in the quasi-stationary approximation. Such a for¬ 
malism is in fact lacking and only partial cases are to be found in the literature 
(see ref. [1,3,11,16]). 

As an application of this general theory we introduce and investigate a 
new design of detectors, normal and superconducting electromagnetic circuits 
which couple to gravitational waves in non-parametric manner. These devices 
are circuit-generalisations of a capacity-device approach invented by Mours and 
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Yvert in 1989. As we shall see in the following onr devices conld reach sensi¬ 
tivities comparable with the sensitivities of mechanical detectors for periodic 
radiation [20,21], 

The paper is organized as follows. In Sect. 2 we discnss, in qnasi- 
stationary approximation, the interaction of electromagnetic helds and cnrrents 
with weak external gravitational helds. In particnlar we discnss nnder which 
conditions the condnctors, i.e. the ions of the lattices of the condnctors, can be 
treated as freely falling in external gravitational wave helds. Sect. 3 presents 
the general theory of normal condncting circnits in the held of external gravita¬ 
tional waves. In this section the parametric natnre of the interaction is shown. 
Sect. 4 rehnes the results of Sect. 3 to superconducting circuits. In Sect. 5 the 
inhuence of dielectric and magnetic matter is investigated. Finally, in Sect. 6, 
we apply our theoretical results to concrete, thermal-noise-dominated detectors 
and derive sensitivity-limits for the detection of gravitational waves. 


2. Electromagnetic Fields and Conductors in Curved Spacetime 

The electromagnetic held equations and the equations of motion for 
charged matter (ions and electrons) in an external gravitational held v = 

0,1, 2, 3) can be derived from the action 


S 


Sem + Sm + S- 


int 5 


( 2 . 1 ) 


where Sem = (l/47rc) J (g = det(( 7 ^,^)) is the action for the free 

electromagnetic held, Sm is the action for the conductors that we will discuss 
below, and Sint = ~ j , with the current density dehned by 

■fi cp dx^ 

^ V-900 dx^ ’ 

is the action for the coupling between the electromagnetic held and the charges 
and currents in the conductor. The electromagnetic held-strength tensor is 
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related with the four-vector potential by — and we have the 

following relation between and the standard electromagnetic held-strength 
entities, Ei and Bi (f = 1, 2, 3), 


F = 


/ ° 

-Fi 

—F2 

—S3 \ 

Fi 

0 

S3 

-S2 

F2 

-F3 

0 

Si 

VFs 

F2 

-Si 

0 / 


( 2 . 2 ) 


By varying the action (2.1) with respect to A^ we get the De Rahm equa¬ 
tions. In the Lorentz gauge, V^A'^ = 0, they read 


- APR}^ = - jp. (2.3) 

P c 

Eq. (2.3) generalizes the inhomogeneous Maxwell equations in the Lorentz 
gauge to curved spacetime; R^i, is the Ricci tensor dehned by 


p _ pa 


where R°‘ is the Riemann tensor 


R^ 




m/3 


+ TZnTl,, 


dxP dxp pp 

and where denote the Christoffel symbols 


pa pM 

^ M/3’ 


r“ = — fif' 

j-iu 


.p f dgpu dgpfj, _ ^ fliy 

\ dxP dxP dxf^ 


(2.4) 


(2.5) 


which enter in the dehnition of the covariant derivative (notations and 
conventions as in ref. [4] pp. 223-4). 

Our aim is to study the interaction of a gravitational wave, emitted for 
instance by a cosmic source, with electromagnetic circuits. In the interaction 
region the metric tensor satishes the Einstein equations in vacuo, R^^, = 0. 
For weak gravitational helds, to which we shall limit ourselves in this paper, 
the metric tensor and its determinant can be written as 

9jii/ F 


.— h 

\/^-l+ 2- 


( 2 . 6 ) 
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where rj^^, = diag(—1,1,1,1) is the Minkowski metric tensor and where 
with << 1, is its small pertnrbation; h = In this ap¬ 

proximation we can £x a particnlar coordinate frame, respectively gange, the 
so-called TT (transverse and traceless) gange [22] specihed by h = 0, /iqm = 0, 
and djh^^ = 0, where the solntions of the homogeneons Einstein eqnations 
take the very simple form: 

h,, = + (2.7) 


where = {ug/c,k) and = 0; A_|_^x are the amplitndes of the two 

polarization states and e_|_ x are the two transverse {eijk^ = 0) and traceless 
(e\ = 0) polarization tensors [22], 

Neglecting second order terms in hgjj one can easily show that the lin¬ 
earized De Rahm eqnations are 


A" = 0. 



(2.9) 


We shall consider onr electromagnetic system in the so-called long-wavelength 
(or qnasi-stationary) approximation in which the freqnency of the electromag¬ 
netic oscillations z/em is mnch smaller than c/d where d is the typical linear 
dimension of the system. This condition is eqnivalent to the assnmption that 
all time derivatives of the potential in (2.9) are negligible with respect to space 
derivatives and also that the total electromagnetic radiation of the system is 
negligibly small (this implies also v/c « 1 where v is the typical velocity of 
the charged particles in the system). If these conditions are satished it has been 
shown ([9] and e.g. [14] §75) that Sem — — (l/2)5'int and thns, the total action 
rednces to 

*5 = ^ / 


5 


( 2 . 10 ) 



In order to get a consistent pictnre of the dynamics of a non-radiating system we 
shall reqnest that the rednced gravitational wavelength Ag/27r is mnch greater 
than d. This entails that the leading term in eq. (2.9) is the one proportional 
to hij. Within these approximations we get 




47r 


- J 

c 




( 2 . 11 ) 


We shall now investigate the matter action Sm- We shall show that if the 
elastic bodies forming the circnits can be treated as free falling in the held of 
the gravitational wave, we can discard the mechanical part of the circnit. It is 
a well known (see [22]) property of the TT system that if freely falling bodies 
are at rest before the arrival of the wave their spatial TT coordinates remain 
constant in time also in the oscillating gravitational held. This fact will greatly 
simplify onr treatment as all the relevant circnit parameter pertnrbations can 
be fonnd in TT coordinates by integration over nnpertnrbed coordinate paths. 

To hnd the free-falling conditions let ns consider the normal modes of a 
solid body as a collection of non-conpled damped harmonic oscillators. This 
problem can be easily solved in the Fermi Normal Coordinates (FNC) frame 
(see [5]). The deformation from the eqnilibrinm shape can be written as a snm 


Sx = Y^ Ar,i(t)'ii)^i{x) ( 2 . 12 ) 

n,l 

The /-label refers to the kind of modes: for instance, as far as cylindrical bars 
are concerned, it refers to the longitndinal, torsional and hexnral modes. The 
n-label refers to the eigenfreqnencies of a given mode. 

For each specihc valne of I the fnndamental mode shall be denoted by 
n = 1. In the following we assnme that only one specihc /-mode is strongly 
interacting with the gravitational wave. For instance, for a spherical body R^i 
is diherent from zero only when 1 = 2. 
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In the presence of a gravitational wave, the coefficients A^i are determined 


by the eqnation (see [4] eqs. (37.42b-c)): 


Anl 4 “ Rnh 

'^nl 


(2.13) 


where 

Rnl = \ j '^nlX^pd^X (2.14) 

(p is the mass density of the body and M its total mass). The solntion in Fonrier 
space (^(w) = ^{t) exp{—iut) dt) is given by 


A-nl (ctl) 





2 

+ iuj/Tnl ’ 


(2.15) 


where is dehned as 



1 

M 


pd^x. 


(2.16) 


For monochromatic gravitational waves of freqnencies Ug » uu we get 


MAt) = i h,j(t) 


(2,17) 


that is to say, in its fnndamental mode the body behaves as a freely moving 
system. This is the operating condition for the interferometric devices nnder 
constrnction in Enrope (VIRGO) and in the United States (LIGO). 

Let ns now consider the opposite case where Ug << uu. We get 

AnAt) = -U‘i,K,{t)(^) (2.18) 

^ \^nl J 

that is the body is practically at rest and the geometrical properties of the 
circnit do not change. 

Finally at mechanical resonance of the fnndamental mode ujg = uju we 
get Aii{u}) = hij{u)Qii where Qu = ttuhTii is the qnality factor of the 
oscillator. This is the operation condition for the resonant antennae already 
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working in Italy (Universities of Rome and National Laboratories of INFN at 
Frascati and Legnaro) and in the United States (Universities of Lonisiana and 
Stanford). 

Now, in order to find free falling conditions in onr case, we examine in 
more details Rni- Comparing eq. (37.45) with eq. (7b) Box 37.4 and taking 
into acconnt dehnition (4) Box 37.4 of [4], it can be reasonably assnmed that 
Rni may have the following behavionr as a f un ction of n 

Rni — -^Rii (2-19) 


This implies that the energy deposited in the n-th normal mode at resonance is 
approximately independent on n (see [4] eq. (37.43)) if we assnme Uni — nuu. 
With these assnmptions the A^i coefficients are given by 


blnz(w) - 


U) 


( 2 . 20 ) 


If we now consider again an incoming gravitational wave with freqnency Ug » 
ujii, different from all the other ujni, then from eq. (2.20) one can see that only 
the fnndamental mode is practically excited. For instance if the hrst harmonic 
excited by the wave is n = 3 (see, e.g. [4], eq. (37.45)) then A^i ~ 
Therefore the body is free falling becanse it follows the incoming gravitational 
wave like a clond of dnst (see eq. (2.17)) and so it is practically at rest in TT 
gange. One can see therefore that, by snitably arranging the fnndamental mode 
freqnency uju with respect to the freqnency chosen for observing gravitational 
radiation, the condition of free falling can be easily fnffilled. In the rest of this 
paper we shall restrict onrselves to physical devices where these conditions are 
satished and therefore we neglect the matter-action term in (2.1). 


3. General Theory of Electromagnetic Circuits 



To establish the behaviour of a circuit in the held of a weak gravitational 
wave we generalize the lagrangian approach of [15]. Let 7 a (a = 1,2,, N) 
be a system of N conducting non-ferromagnetic one dimensional (wires) or 
extended bodies (capacitors). The quasi-stationary approximation allows us to 
neglect surface currents in extended conductors and charged wires. The charge 
over conductor surface and the current Lowing in wires can be treated as not 
depending on spatial coordinates. The relations between, on the one side, 
and the current I and, on the other side, charge density p and the charge dQ 
in an inhnitesimal volume element can be written in any coordinate system as 
(see for instance [14] §90) 


Idx^ 


j' 


(3,1) 


[ dQ = p d^x 

where is the determinant of the three-dimensional metric ^pij = pij — 
goidoj/goo- The equations for / and Q can be found by substituting in the 
action (2.10) the (formal) solution of the De Rahm equation (2.3) for gravita¬ 
tional waves. If we can neglet g^v^a with respect to (i.e. for slowly varying 
gravitational helds) then □ = ~ g^^dadp. In this case we can write 

= “ / K^^[x,x']j''{x)y^-g{x') d'^x’, (3.2) 

where is the retarded Green function, solution of the equation 


Ox Kf,^[x,x'] -47r^^===h^(a:“ - a:'“) (3.3) 

V-gi^) 

In a synchronous reference frame x = {ct,x), where ^gij = gij and —g— ^g, 
we can integrate eq. (3.2) over the time coordinate x'^ and get the generalized 
retarded potential ([4] p. 500) 


Af^{x,t) = -[ K^^[x, x']j’'{x^, a?') a/ ^g{x^, x')d^x'. (3.4) 

c Jv 
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Now we have to introduce the lagrangian £ of a system. ^From (2.10) and (3.4) 
we obtain the total lagrangian of a system of conductors and wires in curved 
spacetime 


C 


1 


' 1 / 


a / '^g{x) d^x / a / x') d^x'x 


'V 


X Kf,^[x,x'] j^{x) j’"(x'^,x') 


(3.5) 


^From the relations (3.1), taking into account that / and Q are only func- 

_Oz 

tions of time in the quasi-stationary-held approximation and that K = 0 in 
synchronous reference frames, we get 


N 


£ _ -^7a 


a,6=1 

1 


dx^ / dx'^Kij[x,x'] 


'lb 


(3.6) 




For quasi-stationary helds charges are distributed over the conductor surfaces 
and so we can replace in (3.6) the charge integrals with surface integrals. 
This leads to 


N 




a,6=1 


2 ^lalbi^)^la ^Ib 


2 ^lalb lb 


(3.7) 


where £ 7 ^ 7 ^ are the generalized coefficients of mutual inductance 


^lalb 


(t) =/ dx^ / dx'^Kij[x,x'], 


'la 


'lb 


(3.8) 


and 


Qla^lalb i^)Q 


lb 




dA 


, dQ 
^^dA 


la 

la 


Koo[x,x'] 


dQib 


(3.9) 


dehne univocally the coefficients of capacity by means of the equations 

X] 6 =l ^lalb^lblc — ^ac- ^From the above formulae one easily recognizes that 
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the action of gravitational fields in synchronons reference frames (snch as TT- 
gange frames) is a change in the geometrical circnit parameters, i.e. the in¬ 
teraction of gravitational radiation with standard electromagnetic circnits is of 
parametric natnre. 

In order to take into acconnt also dissipations of real circnits we go over 
to Ohm’s law written in its covariant form ([8] p. 263); for the cnrrent in wires 
one has 

r - jf3 = a FZ u’^ , (3.10) 


where u'^ is the normalized fonr velocity of the lattice of the condnctor and 
aa is its condnctivity which we assnme constant in the range of freqnencies we 
are interested in. As the lattice has fixed TT coordinates = Jq , the total 
dissipated power in onr circnits by Jonle effect is 

P = [ ^gd^x r 

Jv 

= g^j fix’") fix’"), 

a A. 

which can be rewritten in the form 



N 

P = J2Ry.mil, (3,12) 

a=l 

where Raf) are generalized resistances 


RyAt) = -% [ f^tldx'-, (3,13) 

here tf is the nnit vector ( = 1) tangent to the wire ja and = 

d^x.y^/dl, idR = gijdx'^dxf is the cross-section of the wire. 

In order to get the Enler-Lagrange eqnations from the lagrangian (3.7) we 
have to set np a relation among onr charges and cnrrents. ^From (3.1) and the 
charge conservation relation follows . Hereof the eqnations of motion 
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result in the form 


N r 


E 

b=l 




{t) Qji, -\- 

^lalb {t) Qjf^ + 

^lalb {t) Q 


76 


= 0 . 


(3,14) 


Once we have dehned the dissipation f un ction $ = P/2 we can substitute the 
r.h.s. of the equation (3.14) with d^/dQ-y^ and write (see [13] §1.5) 

N 

^ ^ ^laJbi^)Qlb ^laJbi^)Qlb Qlb ~ ~^lai^)Qla- (3.15) 

6=1 


These equations describe, in TT coordinates, RLC circuits in external gravita¬ 
tional radiation helds when dissipations are taken into account. 

In the weak gravitational held approximation one can easily show (cf. [10]) 

that 

A-,.. = P / V + ^h,a"‘B"\x - x'\ri^, (3,16) 


satishes eq. (2.11) through eq. (3.4), and therefore it gives the right expression 
of the Green function in the linear and quasi-static held approximation 
and for a TT-metric perturbation hy^^it). From eqs. (3.8), (3.9) and (3.13) it 
follows that, in the same approximation, the coefficients and 

read 


^lalbi^) ^ ^lalb ffi^7a76(^) 


= (3,17) 

RiSi) = +’’ 7 .( 1 ) 

where , and ^Rj^ are respectively the usual inductance, capaci¬ 

tance, and resistance coefficients in hat spacetime, dehned as 


^lalb 




dx\ dxb 

la 76 


l^7a ^7bl 


(3,18) 


Up _ X 


^iJla dl 


'la 


dA 


la 


(3,19) 
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Cilo. OA f\jb — 


f r dA^JA'^^ dQ^» 

'la Jib ^Aa l* 7 a ~ *75! 


(3.20) 


Here c^„ 7 b(t) and r^^(t) are the time dependent pertnrbations indnced 

by the gravitational waves on the circnit parameters which can be written as: 

^lalb ~^ij^lalb 


c ^ h • • 

^JaJb '^^J^Ja7b 

T ^ h ' ■ 

I Ja 


(3.21) 


yj ^ _ I I la Jb 

lalb „2 / 1.;^. _ y 

^ Jla Jib I'^la *^76 


dxl dx'^^ 1 

la lb I _ 

- a:', I 2 c 2 


5kid'^d'^\x^^ -x'^^\dxAA^ 


V 2j^J^^dA.„dy 




(3.22) 

(3.23) 


= 

^la 


r 

'la ^la (^^la 


(3.24) 


The coefficients ^Ab^ Aaib ^ 7 a’ defined in (3.22), (3.23) and (3.24), are 
pnre geometric qnantities which play a similar role as the mass-qnadrnpole 
tensor p(x)(x^x-^ — (l/3)S^-^x^Xk) d^x in the mechanical interaction 

between a gravitational wave and a solid body. 

We notice that the eqnations that describe the interaction between grav¬ 
itational waves and electromagnetic circnits are parametric ones. This means 
that gravity is not acting as electromotoric force bnt is changing the circnit pa¬ 
rameters. The effect of the small pertnrbation of the parameters of the circnit 


is a freqnency and amplitnde modnlation (see for instance [18]). 


4. The superconducting circuit 

Below a certain temperatnre T^, called transition temperatnre, some met¬ 
als become snpercondncting. The most impressive properties of this state are 
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the sudden drop of the electrical resistance and the Meissner effect. F. London 
in 1935 (cfr. [12]) suggested a macroscopic theory of the pure superconduct¬ 
ing state which unified these two experimental facts under the same theoretical 
scheme. He supposed that the density current flowing in the superconductor 
was a sum of two currents: the normal and the supercurrent 

j=js+3n- (4-1) 


The two kinds of current are related to the electromagnetic field inside the con¬ 
ductor in different ways. The supercurrent satisfies the two London equations 


curl (Ajg) 






B 

c 

E. 


(4.2) 


Here A is a constant characteristic of the superconductor (its order of magnitude 
is > 10“^^ sec^); B and E are the intensities of the magnetic and electric fields. 
The normal current is connected with the electric field by Ohm’s law 


in = 


(4,3) 


where a is a continuous function of the temperature and has no jump at T = T^. 
The total current j satisfies the Maxwell equations and the relation between 
the two currents turns out to be 


dt aA ^ 

where /? has the dimension of a frequency (/? = 1/aA 10^^ Hz). 
electromagnetic field oscillates at a frequency Uem we have 


I . I ^em I . 

\jn\ = — \3. 


(4,4) 
If the 


(4.5) 


and if Uem << (3 then << |j^|; in this way the sudden drop of the resistance 
is due only to the fact that the normal current (which is the one that dissipates) 
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becomes practically zero. Moreover, from eq. (4.2) and the Maxwell equations 
follows 

^Bi = ^ Bi (4.6) 

where 



is the so called penetration depth. The regular solutions of the eq. (4.6) decrease 
exponentially from the surface of the superconductor to its interior with typical 
length of the order of the penetration depth. Therefore inside large bodies (large 
compared with the penetration depth) we have .6 = 0. 

The London theory of superconductivity is enough for our purposes. In 
fact, as in the preceeding section we didn’t need a microscopical theory of 
conductivity so now we need only the knowledge of the macroscopical behaviour 
of a superconductor that is given by the eqs. (4.1)-(4.3). 

In the following we point out the differences from the ohmic case. The hrst 
step is to write down the equation of an oscillating superconducting RLC circuit 
in flat spacetime. ^From the Maxwell equations we get the energy theorem as 
usual 

div-^(E x B) + ~{B^+ E^) = -HE). (4.8) 

For quasi-stationary conditions the hrst term is negligible. (l/87r)(6^ + E^) is 
the energy density of the electromagnetic held and it has the same expression 
as in the ohmic case. 

The diherences come out when we write explicitly the work done by the 
electric held (see [12] §9), 

The term represents a reversible work (it is in fact the kinetic energy 

density of the supercurrent) while (l/a)j^ is the energy density dissipated by 


15 



the normal current (Joule effect). Another difference arises from the fact that 
if LVern << (3 the Current is different from zero only for a few penetration depths 
below the surface. Therefore we need to be careful when we calculate the kinetic 
energy and the dissipated power. In fact, let us approximate the current as 

(4,10) 

where ro is the radius of the wire (ro >> A); the power lost as Joule heat is 
given then by 

P = f -jldV = RJl (4,11) 

Jv ^ 

where 

/n ^ 27rAro |j^(t)| (4.12) 

and 

Rs = [ (4.13) 

rj is the effective section in which the current Hows. Is takes the value 

rj ^ AtiXtq. (4.14) 


3 (a ^) ~ 3 (t) exp - 


(ro - r) 


In a wire one has Rg ~ (^o/A) R, where R is the usual resistance. For the 
kinetic energy we get 

W,:,^ ^ I £ dV = ^ RJ^ (4,15) 

where for the last equality see eq. (4.4); A is defined analogously to 

Therefore we can write the lagrangian of an RFC superconducting circuit 

as (/ At In-) Q Qs TQn) 


£ + LIJ„ + Rll - ^Ql - - Tq2, (4,16) 
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The dissipation function reads 


(4.17) 

and the relation between In and Ig is given by (see eq. (4.4)) 

is = I31n. (4.18) 

Using the method of Lagrange multipliers (see for instance [13] §2.4) we hnd 

Q + ‘^'JsQtl + ^oQ ~ 0 (4.19) 

where 27 s = Rg/L and cOq = 1/LC. The only difference from the equation of an 
ohmic circuit is in the resistance term: in the superconducting case this term 
is only due to a part of the current (the normal current). 

Making use of relation (4.18) we hnd for the total current a differential 
equation of the third order in time 

T +/3{1 + 2r)i + uli +/3u^I = 0 (4.20) 

where 2r = Rg/f3L << 1. The solution of this equation can be written, up to 
terms of the hrst order in T, as 

I{t) = [B cos (Oq^) + D sin (Oo^)] (4.21) 

where 

Hq =^^ 0(1 ~ r) 

rsi§ (4.22) 

70 = — 

The hrst term is a sort of opening extracurrent which has its origin when, after 
starting the current with an external e.m.f. of frequency we let oscillate 

the circuit freely. The coefficient A is practically diherent from zero only if ooemf 
is very diherent from lvq. In this case the current must change its frequency 
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according to (4.18). The characteristic time of this term, being of the order of 
is too short to be taken into acco un t in a macroscopical theory like this. 
Therefore we can neglect this term and pnt A = 0. In this case it is easy to 
recognize that the following eqnation 

/ + 270 / + = 0 (4.23) 


has the same solntion of eq. (4.20) i.e. an oscillating RLC snpercondncting 
circnit is eqnivalent to a normal one with resonance freqnency Oq and resistance 
Reff=Rsivi//3^. 

For N snpercondncting wires or bodies the electric and magnetic part of 
the lagrangian as well as eqs. (3.8) and (3.9) are the same as in the ohmic case. 
Withont dissipation and snpercondncting kinetic contribntion, the lagrangian 
in cnrved spacetime is still given by eq. (3.7). As far as the dissipation fnnction 
is concerned, recalling that only the normal cnrrents dissipate and that they 
flow only near the snrface, we get 


N 




2 

n~(a 


a=l 


where 


(^) ~ 


9ij 


la 


(4.24) 


(4.25) 


Jla 

in which is the effective cross-section (4.14) of the a-th wire. These eqna- 
tions allow ns to write down the generalization of the kinetic energy 

N 


Wkin{t) — 2 ^ -^S7a ’ 

^ a=l 


(4.26) 


and hnally the lagrangian of a system of N snpercondncting wires or bodies can 
be written as 


AT r 


a,h=l 


1 


1 


2 -^7a76 (^)-^7a -^76 2 ^1 aQ 


1 


N 


76 I 2/1 ^ ■ (4-27) 


a=l 
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The dissipation function is given by (4.24). To get the Euler-Lagrange equations 
from lagrangian the (4.27) and the dissipation function (4.24) we have to know 
the relations between charge and current and between charge and supercurrent. 
For a set of RLC circuits we have and because of the general covariant 

form of eq. (4.2) ( d/dT[K[ji^ — j^UaU^)] = the relation between Qsa 

and Qna is given in the curved case also by 

Q,-,.=PQnj.- (4.28) 


Using the method of Lagrange multipliers we get from (4.27) and (4.24) 


N 


El 

6=1 


^7a76 (^)^76 “i" -^7a76 (^)^7b ^7a76 (^) Q lb 


(4.29) 


and by means of (4.28) we get the equations for the supercurrent 


N 


^ ^ [ ^lalbi^) Q ^(^^lo.lbi^) T ^7^7), (t))(5s76 T Rs 7^ (t)(5s7aT 


6=1 


(4.30) 


S'fb T f^^lalbi^) Qs'fb 


= 0 . 


5. Conductors in the presence of matter 


In this section we generalize our theory to the case in which the capacitor 
is hlled by a dielectric and the inductance has a magnetic kernel. Within the 
usual approximations the matter can be characterized by a dielectric constant e 
and a magnetic permeability n which enter into the constitutive equations (see 
[8] p. 263): 


idctf^FjS T dcti^Sfi T dcxsFp^^ u 


d (.dctP^iS T 9ct~fH§p 9asH 


u 


a 


(5.1) 


where is the macroscopic 4-velocity of the medium. With these dehnitions 
the Maxwell equations become 


F F F = 0 

^ in',<7 I ^ l/(T,ld I ^ <711,ly ^ 

4:71 

= —r. 

’ c 


( 6 . 2 ) 
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Under free-falling conditions = Sq, eqs. (5.1) rednce to 


Hok — e Ff 


Ok 


Fik H, 


ik • 


(5.3) 


In this way the second pair of Maxwell eqnations becomes 

.•0 


47r / 

e c 

drr .u 
F% = — 

’ c 


(6.4) 


The action of the free electromagnetic held reads S, 
(l/47rc) / 

With these notations eqs. (3.18) and (3.22) become (see [15] 

Or d 


dxl, dx'4, 


>r — 

~ „2 I I "*.7 \ 

^ Jla Jib -^76 1 


30-33) 

(5.5) 


[ _22_ Ih _ — [ [ Skid d'^lx -X \dx^ dx 

'7a 7b ^2 / / I™ _ 7 ./ I ^ 9r>2 / / I -^761 7a “-^76’ 

27,, J7(, |•^ 7 a -^76 1 27,, J7b 


\ jj _ d_ 

r> 


7b I 7a 7b’ 

(5.6) 


while (3.20) and (3.23) become 

/nia OA ^7b _ 
^lalb^ 


dQ^“ dA^JA' dQ'^‘^ 


_ _76_ 

^ Aa Jib ^^la l*7a ~ *7 I *^^7 


Q'^^FQ'"- = 4: 


2e 


d(5^“ dQ^'’ Q/iQ/7i 

-tt —d a Ax ^ 

dA^ dAL 


*76 I *^^ 7 a 


76 


(5.7) 


(5.8) 


'la -Jib '^"^7a ™*7b 

Therefore the presence of matter with scalar dielectric or magnetic properties 
changes only the dehnition of the system parameters; they are increased by e 
and n factors. This is important becanse, in order to obtain great capacitances 
one can nse snitable dielectrics: standard commercial types of capacitors can 
reach capacities of abont 10 mF. 


6. Application and Discussion 

Let ns now consider the device of Fignre 1, located in the x — y plane 
with X = x^ and y = x‘^, and evalnate its response to a periodic gravita¬ 
tional wave nsing realistic parameters. This circnit can be described as a set 
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of six conductors: two wires labelled 71 (the inductance and resistance) and 72 
(idealized connection wire) which connect four plane conductors labelled with 
73 , 74 (to form the capacitor Ci) and 75 , 75 (the capacitor C 2 ) with charge 
Qx Q Q'y 4 Qi Q Q Q 2 Q and Q^q Q 2 Qi and 
Q 2 are the constant electrostatic charges on the plates of the two capacitors in 
the absence of gravitational waves when there is no current flow which satisfy 

Qi + Q 2 =Qo, 


( 6 , 1 ) 


Qi/Ci — Q2IC2 — 0 . 

The connection between the capacity of a condenser and the capacitance coef- 
hcients of the conductors i and j reads (see [15]) 

C-^ =Cu-2C,j+Cjj. ( 6 . 2 ) 


The lagrangian of this system can be written as 


, HQ-Qir HQ+Q2r 

2 ^ 2 C'i(t) 2 C 2 it) ■ 

Taking into account dissipation, the equation of the circuit becomes 

Qi Q 2 


L{t)Q + L{t)Q + R{t)Q + 


1 1 
+ 


Q = 


(6.3) 


(6.4) 


Ci(t) ' C2(t)J~‘ Ci(f) (72(f)' 

Let us now consider the case of a weakly coupled gravitational-wave held. Set¬ 
ting 


_ . ^ ^33 

^ + ,x — (^ + ,x)ij 

^+,.=-ie+,.h{x\\-‘^X% + xl^2)C 

d+,X — + 


(6.5) 


we can separate the polarization and angular dependences of the interaction 
(hgure pattern). In this way we obtain 

L{t) = L[1 + eLit)] = L[1 + a+h+{t) -S axhxit)], 

( 6 . 6 ) 


= ^[l-ec(t)] = ^[l-K+h+{t) -^^x/^x(f)], 


R{t) = R[1 + eR{t)] = R[1 + p+h+{t) +pxhxit)]. 
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For the sake of simplicity let us put 

Cl = aC 


Uq 


C2=C 
a + 1 1 
a LC' 


With this dehnition of the factor a we obtain 


Qi I 1 Qo 

a + 1 


Q2 — — r~rQo 

a + 1 


and the equation of the circuit can be written as 

(1 + €l(«))<? + 27 (1 + ^ + €«(«)) <? + (1 - j Q = 

2^ a(eci(f)-ec2(f)) 

- ‘"oQo--■ 


As the time dependent coefficients in the l.h.s. are very small compared with 
unity their unique effect will be an amplitude and frequency modulation of the 
unperturbed charge motion 


Q + 2'yQ + UqQ = v+h+{t) + (t) 


( 6 . 10 ) 


where 


'^+,x — ^oQo 


2 „ — ^ 2 ,+,x) 


( 6 . 11 ) 


-T.x (a+ 1)2 ' ' 

In equation (6.10) the interaction between the gravitational wave and the circuit 
is non-parametric. The dropped modulations are so small that they cannot be 
observed being a second order effect with respect to the solution of (6.10). If we 
consider a periodic gravitational wave h+^x{t) = A_|_^x cos {ugt + (f)+^x)j then 
the solution of (6.10) can be written as 


Q{t) — (5_|_(f) + Qx (f) 


( 6 . 12 ) 


Q+,x (f) — 


2A+.,Q„Q ^ + 


(a + 1)^ 


2_,,,2'>2 Un 


K--?) 
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where Q = c<;o/47 is the quality factor of the circuit and 

h = arctan —(6.13) 

In this system the effect of the gravitational wave is to cause a current I{t) to 
flow in the circuit with 


2AuoQoQ + , a , 

- 1 . „ cos [Ugt + 0 + 6). 


{a + ly 


^ _j_ '• 0 g> 


(6.14) 


47^ 


'T71T- 


Near resonance (i.e. (wq — << (^7 )^) we have 

r a(^^l-^^2) 27lwo<5oQ , J. , r ^ 

I near res -i \ 9, / \ 2 COS [yjgZ (p Onear res J 


(a + 1 )^ 


l + 2 (^^ 


(6.15) 


where hnear res — . This device behaves formally like a Weber bar with 

the initial static charge Qo corresponding to the rest length of the bar. The 
device has two advantages: hrstly, it can be easily tuned to receive radiation 
from periodic sources like pulsars and binary stars (i9g > few Hz); and secondly, 
the measurement is a null measurement (the expected current versus zero). 

Let us now consider the superconducting case. Using the same notations 
eq. (6.4) becomes 


L{t)Q + L{t)Q + Rs{t)Qri + 




Q 


Qi 


Q 2 


(6.16) 


For weak gravitational waves we have 


(1 + ^L{t))Q + ^L{t)Q + 27 s [1 + ^R{t)]Qn + 1 

^oQo ' 


,2^ a{ecyt) - ecyt)) 


(a + 1)^ 


eci (t) + aec2 jt) \ ^ 
a + 1 J 


(6.17) 

Again the time dependent coefficients in the l.h.s. of equation (6.17) produce 
only a very small frequency and amplitude modulation of the solution of the 
following equation 


Q + 2'ysQn = v+h+{t) +Vxhx{t) (6.18) 
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where 


,,+.x = SJ ?(1 + 2r)Qo^hkti2L_pi£l, (6,19) 

In equation (6.18) the gravitational-wave electromagnetic-circuit interaction is 
purely non-parametric. If we use relation (4.28) in (6.18) we hud the equation 
of the supercurrent as 


Q g (3 {1 + 2r)(5s + 3- 2r)(5s + 3- 2r)(5s — P[v^hj^{t) + Vxhx{t)]- 

( 6 . 20 ) 

If h+^x{t) = Aj^^x cos {ugt + (f>+,x)j then the solution of (6.20) can be written 
as 


Qsit) — Qs,+it) + < 3 s,x (f) 


( 6 . 21 ) 


2-4+,x(l + 2r)Q„Q.„ no + 


where 


and 


(a+ 1)2 T{ug 


(i + 2 r)( 02 _(^ 2 ) 

t<3,n 0 — 


OJn 


270^^9 (^1 + 


27od 


T{ujg) — 1 + 


02 - 


u 


Q - ^0 _ 
2 e// - XT - 


270/3 

/3 


- I a/i + tan2 5 


(3^ L 


( 6 . 22 ) 


(6.23) 


(6.24) 


470 4 rOo 2O0 Rs 
is the quality factor of the superconducting circuit which is several orders of 
magnitude greater than that of the ohmic circuit. In this way the gravitational 
wave gives rise to a supercurrent 


<^(^1 ~ ^2) 2 A(i + 2 r)Oo( 5 oQe// / + 1 7 1 

I sit) = -^-2+ COS [Ugt + (j) +5) 


(a + 1 )" 


T{ujg 


(6.25) 


in which, for the sake of simplicity, we have considered only one polarization 
state and dropped the indices +, x. Near resonance (i. e. (Oq — C. 2)2 « ( 2 yo) 2 ) 
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we have 


‘■s near res 


, . a{ni — K2) ‘^AQ.oQoQeff / . , / , X 

(t) = - 77 - -T2 cos {ugt + (j) + 5, 


near res > 


(6.26) 


( g . + 1) 1 2 ^ ^0 ^ 

where hnear res — ■ To get the total current we have to use (4.28) for the 

normal current, 


In(t) = “i"*, 2r)»00,Qe// ^ j 

/t (a + 1 )^ ^"(‘^g) 

We see that this current (in agreement with (4.4) and (4.5)) is smaller by a 
factor oJg/jS and it is 7 r /2 out of phase with respect to the supercurrent. The 
response of a superconducting device differs from the response in the ohmic case 
because of different quality factors. 

Let us now calculate explicitly the parameters of the circuit. The eqs. 
(3.22-24) and (3.18-20) or (5.5-6) and (5.7-8) give: 




(6.28) 


We note that 


d' d' -x'^,\ = 


[y-y'Y + {z-z ') 
l*7a ~ ^76 


/^2 


etc. 


(6.29) 


so that 




dQ-<‘ dQ-'^ (x' - I'M 


n \2 

dA.., dA' 


/ dA dA' W -t' 

7a Jib '^-^la, “^76 I•^7a *^76 I 

(6.30) 


(no summation over i). 

Let us now consider the case of a capacitor made of two circular plates set 
perpendicular to the x axis and at a distance d apart. If 7 ^ = 75 then 


X. 


11 

7a 7a 


= -a 


7a 7a 


(6.31) 
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and 


X 


22 

7a 7a 


= X 


33 

7a 7a 


-^(7 

~ ^'^7a7a- 


(6.32) 


If 'ja ^ lb, assuming that d is much smaller than the diameter of the circular 
plates, by means of (6.28) and the symmetry of the problem, we approximately 
get 


^7a7b 


_ n 

^lalb 2 


X 


22 

7a 76 


X 


33 

7a 76 


—C 


(6.33) 


which is valid also for 7^ = 75 as can be seen by comparing (6.33) with (6.31-32). 
Similarly one can show that X7^75 = 0- 

Now let us assume that the wave comes from the z — axis perpendicular 
to the plane containing our circuit (see Figure 1). In this case we have for the 
polarization tensors 


(e+)^!y 



0 

0 


0 

1 

0 

0 

0 

0 

-1 

0 

Vo 

0 

0 

0/ 

/o 

0 

0 



( \ _ 0 0 10 
(ex)/.^ -'0100 

Vo 0 0 0 / 

With reference to eq. (6.5) let us calculate one hnds 


(6.34) 


(6.35) 


^+ = -(e+).,(xl^-2xl^ + X*Oc = 

^ ^ (6.36) 

= [(x?? - x}9 - 2 (x?^ - xia + (xi - x^y] a 

Taking the x axis parallel to the plates of Ci and the y axis parallel to the 
plates of C 2 we obtain 


( k+)i = 1 , (^^^)^=0 

(«+)2 = - l ’ (« x ) 2=0 

Substituting (6.37) in (6.15) and (6.26), and taking a = 1 one gets 

Ires = A^uoQo Q cos {uot + (f )), 

Ires.s = A+floQoQeff COS {uot + (j)). 


(6.37) 


(6.38) 
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We now estimate the minimal gravitational amplitnde which can be detected 
by onr device. Once the device is isolated from the ontside world by means 
for instance of mechanical hlters, as those nsed in gravitational bar detectors 
and Faraday cages, the main noise sonrce is the thermal one, as well for the 
mechanical part of the detector as for the electric part, e.g. the back-action 
noise of the ampliher, at cryogenic temperatnre and for the best amplihers, is 
at most of the order of magnitnde of the thermal noise. In the following we 
shall derive formnlae for the thermal noise in the electric part of the detector. 
The thermal spectral noise in the mechanical part {261/1 = h, see eqnation 
(2.17)) is given by ^JlQkBT/MrPuj^ where T, M, r, and I denote temperatnre, 
effective mass, life-time of the fnndamental mode, and effective length of the 
detector, ks is the Boltzmann constant and u denotes the angnlar freqnency of 
the gravitational wave, e.g. see [24]. Assnming the nnmbers T = 4 K, M = 10^ 
kg, r = 10 sec, / = 10 m, and u/27r = 60 Hz yields abont 2.1 x 10“^^ for h, 
measnred throngh a period of 4 months. This nnmber is as small as the smallest 
nnmber in the eqnation (6.45) below, i.e. for the snpercondncting circnits in 
qnestion the mechanical part has to be adjnsted to those nnmbers. 

We now make the hypothesis that the only resistance in onr circnit is that 
of onr scheme in Fignre 1. The mean sqnare flnctnation of the voltage noise is 
therefore given for tabs » Q/^o by ([25] p.288) 

< >= ^ (S.39) 

be// 

where Ceff is the total capacitance of the capacitors. The mean noise energy 
is 

Cef, < V/ >= keT (6,40) 

The energy dissipated by the ohmic and snpercondncting cnrrent prodnced by 
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the gravitational wave is 


W = j RPdt = -Rl^tobs 
Ws= I Rslldt = ]^RsllQtobs 


(6.41) 


where tobs is the observational time and /q the amplitnde of the cnrrent /. 
Eqnating (6.40) and (6.41) one hnds the minimnm detectable cnrrent after an 
observation time tabs, 

Io = 


2/cfj T 


InO — 


R tobs 

2kB T 


(6.42) 


Rs tobs 

At resonance this implies (by means of eqs. (6.7), (6.15) and Q = uJoL/{2R) 
for the ohmic circnit and (4.22), (6.24) and (6.27) for the snpercondncting one) 
that 


A, 


1 


(a + 1)^ 


keT 


Vo y - K2y QCuotobs 


A, 


(a + 1)^ 


keT 


(6.43) 


noise s 


Vo V - K2Y QeffCuJotobs 

where Vo = Qo/C . If a = = —K 2 = 1 then 


A, 


A. 


2 

2 


keT 


QCuotobs 


keT 


(6.44) 


Qef fCuJotobs 


Taking Vq = 10® Volt, T = 4 K, (7 = 10 ^ F, wq = 27r x 60 rad/sec, tobs = 4 
months = 10^ sec, Q = 10^ and Qeff = 10® one has 


A„„,sp — 7.6 X 10 


-22 


A, 


noise s 


~ 2.4 X 10 “^^ 


(6.45) 


Finally we compnte the order of magnitnde of the normal and snpercondncting 
cnrrent, prodnced by a gravitational wave with amplitndes A_|_ = 2 x 10 and 
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A+ = 5 X 10 which correspond respectively to a signal to noise ratio eqnal 


to 2; we hnd 

/n ~ 7.5 X 10“^^ A 
Is^ 2x 10“^^ A 


(6.46) 


that are within the possibilities of actnal measnrements with SQUIDs. If we go 
down with the temperatnre in the nltracryogenic zone (~ 40 mK) we obtain 

~ 7.6 X 10“23 / ~ 7.5 X 10“^^ A 


(6.47) 

A^o^se s ^2.4 X 10-24 ~ 2 X 10-42 A 

Fnrther improvements can come either from different arrangements of the circnit 
elements (for instance one can consider a series of many capacitors). As far as 
the last qnestion is concerned, today’s limitations in the snpercondncting case 
are dne to the energy dissipation in the material which acts as a monnting of 
the circnit. For this reason we have set Qeff = 10® bnt there are no theoretical 
limitations in thinking of qnality factors higher by several orders of magnitnde: 
it is in fact only a technological problem. We think that a serions research in 
this held wonld be very nsefnl. 

It is to be expected that when the device is operating near the mechanical 
resonance the sensitivity shonld increase. However this sitnation deserves a 
more detailed investigation. 

Recall now that we have calcnlated the cnrrents in TT-gange. In ac¬ 
tnal experiments we measnre cnrrents in the laboratory reference frame (FNC). 
Becanse of the charge conservation we can write 


dy^ = dx^ 


(6.48) 


where are the FNC coordinates while are the TT coordinates. As x^ — 
— \h-’^y'^y^ we have 


dx^ 


1 - 


- ^h\fy^dy^ 


(6.49) 
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where the superscript (i) is the i-th derivative of h with respect to its argument. 
Therefore the relation between the current in the two gauges reads 


jFNC ^ jTT I ^ 


i 'i jy I 'l^dyy 


(6.50) 


In our approximation both hrst order corrections are negligible and so 


jFNC ^ jTT 


(6.51) 


We point out that our theoretical limits in the magnitude of the 
gravitational-wave amplitude detectable from periodic sources like the Crab 
pulsar are better than those experimentally obtained by a Japanese research 
group (A < 5 X 10“^^ [20]; see also [21]). Therefore, a superconducting device 
at liquid helium temperature could give new limits for the emission of gravita¬ 
tional waves from this pulsar. As to vary an electromagnetic frequency is much 
simpler than varying a mechanical one, our detectors can be easily adjusted to 
any gravitational-wave frequency measurable on Earth. 
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Figures 

Fig.l — A scheme of the RLC circnit describrd in the text which can be 
nsed as a detector of gravitational waves. A static charge is distribnted on the 
plates of the condensers Ci and C 2 - 
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